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Abstract 

A number of different fits to solar neutrino mixing and mass square difference were 
performed using 1496 days of Super-Kamiokande-I's solar neutrino data. These data 
select two allowed areas at large neutrino mixing when combined with either the 
solar 8 B flux prediction of the standard solar model or the SNO interaction rate 
measurements. A global fit combining SK data with the solar neutrino interaction 
rates measured by Homestake, SNO, Gallex/GNO and SAGE prefers a single allowed 
area, the Large Mixing Angle solution, at the 98.9% confidence level. The mass 
square difference Am 2 between the two mass eigenstates ranges from about 3 to 
19 x 10 _5 eV 2 , while the mixing angle 9 is in the range of tan 2 9 ~0. 25-0. 65. 



1 Solar Neutrino Oscillations 



While a number of previous solar neutrino 
experiments [1-4] have measured a smaller 
neutrino interaction rate than that predicted 
by the Standard Solar Model (SSM) [5], in 
1996 Super-Kamiokande (SK) began collect- 
ing data of unprecedented precision and quan- 
tity. In addition to a reduced interaction rate 
of solar 8 B neutrinos (originating from (3 + de- 
cay of 8 B nuclei) with electrons in the detec- 
tor, SK measured with high precision the spec- 
trum of the recoiling electrons as well as time 
variations of this reduced rate [6]. 

Neutrino flavor oscillations are able to explain 
all observed solar neutrino interaction rates, 
including the rate difference between SK and 
SNO [7]. The large neutrino flavor mixing be- 
tween the second and third generation inferred 
from atmospheric neutrino data [8] in con- 
junction with the absence of an oscillation sig- 
nal in the CHOOZ reactor neutrino experi- 
ment [9] (which is sensitive to oscillations in 
a similar region of mass square difference) re- 
quires the neutrino flavor mixing between the 
first and third generation to be small. Solar 
neutrino oscillations can therefore be approx- 



imated by a two-neutrino description with the 
parameters (mixing angle) and Am 2 (dif- 
ference in mass square). The mass eigenstate 
nil { m 2) is chosen to have the smaller (bigger) 
mass Q. Therefore, Am 2 = m\ — m\ is always 
positive and the mixing angle ranges between 
and 7r/2. For Am 2 between 10 _8 eV 2 and 
10 _3 eV 2 , matter densities in the sun and the 
earth can strongly affect the oscillation proba- 
bility. Consequently, the symmetry of the vac- 
uum oscillation probability around 9 = tt/4 is 
broken. For 9 < 7r/4, resonant enhancement of 
the oscillations (MSW effect [10]) can occur in 
the sun, while for 9 > tt/4, an anti-resonance 
can suppress the oscillations. Since matter- 
dominated solutions usually don't occur for 
9 > 7r/4, we have not shown this "dark side" 
of the parameter space in previous reports. 
Below 10 _9 eV 2 , the oscillation probability is 
more affected by the oscillation phase than by 
matter effects. A logarithmic scale of the vari- 
able tan 2 9 illustrates the mirror symmetry of 



Alternatively, it is possible to call the mass 
eigenstate mi (1712) which has the larger cou- 
pling to electrons (muons). In that case, 9 ranges 
between and 7r/4 and Am 2 might be nega- 
tive. The oscillation probability remains invari- 
ant under the transformation 9' = ^ — 9 and 
Am 2 ' = -Am 2 . 
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Fig. 1. (a) Allowed areas from other experiments measuring the solar v e flux. The shaded area uses 
only Homestake data, while the hatched area uses only the SNO charged-current rate. Overlaid (inside 
dashed lines) is the region allowed by Gallex/GNO and SAGE, (b) Allowed regions from a combined 
fit to these charged-current rates. All contours in this and other figures are 95% C.L. 



the vacuum region (below ~ 4 • 10~ 10 eV 2 ) at 
the maximal mixing line. Figure 1 (a) shows 
the experimentally allowed areas using vari- 
ous measurements. 

In addition to the suppression of the v e rates, 
neutrino oscillations also produce other neu- 
trino flavors. Furthermore, they can induce 
distortions of the neutrino spectrum and time 
variations of the solar neutrino fluxes. In the 
MSW region, the time variations arise from 
matter effects inside the earth (daily varia- 
tions). In the vacuum regions, the time vari- 
ations are the consequence of the change of 
the oscillation phase due to the yearly varia- 
tion of the oscillation baseline - the distance 
between the sun and the earth. 

In this paper we report on the new SK data 
set of 1496 days, which represents the comple- 
tion of the analysis of all the data from Super- 
Kamiokande-I. Previous SK energy spec- 
trum and time-variation measurements ruled 
out large regions of the oscillation parame- 
ter space allowed by other experiments [11]. 
When combined with the SSM, these mea- 
surements also restricted tan 2 6 to be bigger 



than ^0.2 disfavored oscillation to a sterile 
neutrino. Our final data set confirms our pre- 
vious measurements and further constrains 
the range of possible oscillation parameters. 

Even stronger constraints can be obtained by 
combining SK results with the rate measure- 
ments of other experiments. In particular, the 
combination of SNO's measurement of the so- 
lar 8 B v e interaction rate [12] with the SK rate 
yields additional information [7,13] about 8 B 
flavor conversion, since SK also has sensitiv- 
ity to other active neutrino flavors. If the v e 
survival probability is moderately small for 
8 B neutrinos (less than about 50%) then SK 
can observe a significant excess due to the 
presence of other active species in its precise 
measurement of the elastic scattering rate. A 
similar signature can be observed in the com- 
bination of SK and the Homestake rate. In 
this case, however, the lower energy thresh- 
old of Homestake and its resulting sensitivity 
to other than 8 B neutrinos lead to an ambi- 
guity between spectral distortion and appear- 
ance signatures and therefore to a larger un- 
certainty in the observed excess. 
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Fig. 2. (a) SSM flux independent excluded areas (gray) using the SK zenith spectrum shape alone 
overlaid with the allowed regions of Figure 1 (b) in light gray. The overlap of both is shaded dark gray, 
(b) Allowed areas using only SK data and the SSM 8 B neutrino flux prediction. Both allowed regions 
indicate large neutrino mixing. 

2 Parameter Estimation for Two- All regions have similar x 2 (though the LOW 
Neutrino Oscillations fits are slightly worse with x 2 > xt in + 3.1). 



Without SK data, the strongest constraint 
of the parameters is a x 2 fit to the charged- 
current neutrino rates from Gallex/GNO, 
SAGE (combined into a single "Gallium" 
rate), Homestake ("Chlorine") and SNO ex- 
periments [||]. The light gray areas of Figure 1 
(b) use only these three rates and their SSM- 
based predictions. The regions can be under- 
stood as the "overlap" of the three allowed re- 
gions of Figure 1 (a). Above Am 2 ps 10~ 5 eV 2 
near tan 2 9 = 1 (maximal mixing) all three 
regions meet and form the large mixing an- 
gle (LMA) solution. The small mixing angle 
solution (SMA) is the "crossing point" of the 
Gallium region and the Chlorine region and 
lies between Am 2 ps* 10 _6 eV 2 and lCT 5 eV 2 at 
tan 2 9 ps 10~ 3 . The LOW solution is the large 
Am 2 part of the extensive region(s) between 
m lCT 7 eV 2 and 10~ 9 eV 2 , while the lower 
part of this region is called the quasi-vacuum 
(quasi- VAC) solution. Finally, the vacuum 
solutions (VAC) are found below ps !CT 10 eV 2 . 



SK data breaks the degeneracy in x 2 of these 
regions by analyzing the shape of the "zenith 
angle spectrum" [11] which combines spec- 
trum and daily variation analyses. Here we 
utilize the zenith angle spectrum from the en- 
tire SK-I data set collected between May 31st, 
1996 and July 15th, 2001 (1496 live days). The 
numerical results are listed in Tables 1 and 2. 
There is no indication of spectral distortion, 
the x 2 to an undistorted spectrum is 4.6 for 
7 degrees of freedom (71% C.L.). No signif- 
icant daily variation is found, the day/night 
rate asymmetry is A DN = a5 ° D a ^ ht) = 
-0.021±0.020(stat)t£;^(syst). The gray ex- 
cluded regions of Figure 2 (a) are the result of 
the zenith angle spectrum fit. The SMA and 
VAC regions are excluded since they require 
a distorted 8 B neutrino spectrum. The lower 
part of the LMA and the upper part of the 
LOW solution predict daily variations and 
are therefore disfavored. The excluded streets 
are independent of the SSM neutrino fluxes. 
While the SK zenith angle spectrum shape is 
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Table 1 

SK rate and uncertainty for eight energy bins. The rates and statistical and systematic uncertainties (of 
the spectrum shape) in the third column are presented in units of SSM expectation. These systematic 
uncertainties are assumed to be uncorrelated in energy. The uncertainties in the fourth (uncertainty 
of the 8 B neutrino spectrum), fifth (uncertainty of the energy scale of 0.64%) and sixth column 
(uncertainty of the energy resolution of 2.5%) are fully correlated in energy (but uncorrelated with 
each other). The combined uncertainty (last row) is based on the sum of all bins of the zenith angle 
spectrum. The combined rate has an additional systematic uncertainty of usk = -2 6 % ffl (excluding 
8 B shape, energy scale and resolution), which was added to the uncorrelated systematic uncertainties 
of the spectrum shape. 



Bin 


Range [MeV] 


Rate±stat±syst [SSM] 


B Spectrum 


E-Scale 


E-Resol. 


1 


5.0-5.5 


0.467±0.040+qq}J 


+0.04% 
-0.02 70 


+0.09% 
-0.01 /o 


+0.23% 
-0.21 % 


2 


5.5-6.5 


0.458±0.014±0.007 


±0.1% 


±0.2% 


±0.2% 


3 


6.5-8.0 


0.473±0.008±0.006 


±0.4% 


±0.6% 


±0.2% 


4 


8.0-9.5 


0.460±0.009±0.006 


±0.9% 


±1.3% 


±0.1% 


5 


9.5-11.5 


0.463±0.010±0.006 


+1.7 % 
-1.6 % 


+ 2.5 (y 

-2.4 % 


±0.2% 


6 


11.5-13.5 


0.462±0.017±0.006 


+3.1% 
-2.7 /o 


+4.4% 
-4.1 /o 


±1.1% 


7 

8 


13.5-16.0 
16.0-20.0 


0.567±0.039±0.008 
0.555±0.146±0.008 


+5.1% 
-4.2 /o 
+7.7% 
-5.6 % 


+7.0% 
-6.4 % 
+10.6% 
- 9.6 /o 


±3.2% 

+8.4 <v 
-7.9 /o 


Comb. 5.0-20.0 


0.465±0.005 Iqq}^ 


+1.2% 
-1.0 % 


+1.7% 
-1.6 % 


±0.3% 



Table 2 

Subdivision of bins 2-7 according to the solar zenith angle 9 Z . The range of cos# 2 is given for each 
bin: cos# 2 < is 'Day' and cos# 2 > is 'Night' ('Mantle' and 'Core'). The rates are given in units of 
0.001 x SSM. Only statistical uncertainties are quoted. All systematic uncertainties (see Table 1) are 
assumed to be fully correlated in zenith angle. 





Day 






Mantle 






Core 


Bin 


-0.974-0 


0.00-0.16 


0.16-0.33 


0.33-0.50 


0.50-0.67 


0.67-0.84 


0.84-0.974 


2 


453±20 


442±53 


379±49 


472±45 


522±45 


503±49 


426±52 


3 


474±12 


530±34 


506±30 


438±26 


478±26 


451±28 


439±31 


4 


448±13 


463±36 


470±33 


462±29 


509±29 


461±32 


451±35 


5 


453±15 


449±40 


502±38 


451±32 


473±32 


477±35 


483±40 


6 


477±25 


509±67 


351±55 


391±49 


498±53 


434±56 


521±64 


7 


511±54 


570±150 


831±167 


694±131 


665±127 


441±118 


469±131 


1-8 


459.9±6.7 


483± 18 


476±17 


451±15 


496± 15 


467±16 


456±17 
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Fig. 3. Solar neutrino oscillation parameter constraints from combined fits using the SK zenith spectrum 
and rate and other experimental rates, (a) SK and SNO combined fit. The dashed-dotted lines indicate 
the allowed area if the SNO neutral-current rate is removed from the fit. (b) The combined fit to SK 
data and all rates favor the LMA solution. The best fit point is indicated by the asterisk. All fits are 
independent of the SSM 8 B flux and hep flux. 



best described by quasi- VAC solutions, LMA 
fits have similar x 2 . 

Since there is no indication of a distortion of 
the SK zenith spectrum, an analysis of the 
zenith spectrum shape alone can only yield 
excluded regions. However, if the 8 B flux is 
known, the SK rate can be predicted. Then, 
the zenith spectrum shape can be combined 
with the SK rate into a single measurement. 
Figure 2(b) shows the allowed regions using 
SK rate and zenith angle spectrum in combi- 
nation with the 8 B flux prediction and uncer- 
tainty of the SSM. Unlike other experiments, 
whose allowed areas are shown in Figure 1 (a), 
SK allows only large mixing at 95% C.L, and 
there are two allowed regions: the LMA and 
quasi- VAC solution. 



3 Combination of Super-Kamiokande 
Data with Other Experiments 

A comparison between the SK rate and 
SNO's charged-current rate [12] yields an- 



other 8 B flux constraint which is indepen- 
dent of the SSM, so a combination of the 
SK rate and zenith angle spectrum with the 
SNO charged-current rate need not rely on 
any neutrino flux prediction [13]. The SK 
rate of 0.465tg:8i3 x SSM (see Table 1) ex- 
ceeds the SNO charged-current rate [12] of 
0.349±8;g|i x SSM, by about 4.5(7. If this is in- 
terpreted as an appearance signal of other ac- 
tive neutrino flavors, these flavors contribute 
about 25% to the SK rate and 70% to the 
8 B flux (the v — e elastic scattering cross sec- 
tion for the other flavors is six to seven times 
smaller than for v e ). As a consequence of this 
appearance signal, the two allowed regions 
(dashed-dotted lines) of Figure 3 (a) result 
from a combined fit: a LMA and a quasi- VAC 
solution. If the SNO neutral-current rate is 
used as well, the quasi- VAC solutions are dis- 
favored as seen in Figure 3 (a) (gray shaded 
areas) . 

To constrain the quasi- VAC region below the 
LOW solution, it is necessary to add the Gal- 
lium [3,4] (74.8±|;5 SNU or 0.584l°;°|°xSSM) 
and Chlorine rates [1] (2.56 ± 0.23 SNU or 
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0.337 ± 0.030 xSSM) to the fit. Such a global 
fit is still not influenced by the 8 B and hep 
neutrino (arising from 3 He+p^ 4 He+e + +z/ e ) 
flux predictions of the SSM, which suffer from 
the largest uncertainties. However, it relies on 
the other SSM neutrino fluxes, in particular 
the 7 Be flux (10% uncertainty) and the neu- 
trino fluxes of the CNO cycle (~ 20% uncer- 
tainty). Those fluxes contribute [5] about 15% 
( 7 Be) and 6% (CNO) to the Chlorine rate and 
27% ( 7 Be) and 7% (CNO) to the Gallium rate 
in the SSM. The allowed areas of this global 
fit shown in Figure 3 (b) looks quite different 
when compared with the fit without SK data 
displayed in Figure 1 (b): only the upper part 
of the LMA survives. 



Table 3 compares the four smallest local 
minima of the x 2 describing the global fit. 
The best fit is located in the upper LMA 
The 8 B flux resulting from this fit is 
somewhat higher than expected by the SSM 
(5.05±o:il x 10 6 /cm 2 -s) but well within the un- 
certainty. The hep flux is considerably higher 
than expected by the SSM (9.3 x 10 3 /cm 2 -s); 
however, the uncertainty of this prediction is 
thought to be very large. The fit agrees with 
the SK zenith angle spectrum moderately well 
and easily accommodates the Gallium, SK 
and SNO rates, though the predicted Chlo- 
rine rate is about 2a too high. This worsens 
the otherwise very good best fit x 2 ■ 

The quasi- VAC solution describes the SK 
zenith angle spectrum slightly better than the 
LMA does, however, the rates do not fit well. 
The resulting best fit 8 B flux falls 1.5<r short 
of the SSM prediction. The surviving (at 99% 
C.L.) LOW solution fits the rates about as 
poorly as the quasi- VAC solution. The SK 
zenith angle spectrum fits LOW worse than 
it does either LMA or quasi- VAC. The lack of 
zenith angle variation in the SK data reduces 
the Am 2 (usually around 10~ 7 eV 2 ) and wors- 
ens the LOW best fit which is already under 
pressure from the rates. 




0000000000 

10 20 30 40 50 60 70 80 90 

e 

Fig. 4. Magnified view of the LMA solution. 
Note the linear scale of the x-axis. The 95% C.L. 
contours are obtained from the global fit with 
(dashed, dotted and dashed-dotted lines) and 
without (gray shaded area) using the SSM 8 B flux 
calculation. The dashed-dotted lines include all 
experimental information, while the dashed (no 
SNO rates) and dotted lines (no radio-chemical 
rate measurements) do not. 

At around the 3.5<r level a small mixing angle 
solution appears. The mixing angle is about 
a factor of four smaller than the usual SMA. 
For any given 8 B flux, such a small mixing 
angle region can be defined by the crossing 
of the Gallium and the Chlorine allowed ar- 
eas; those two rates therefore automatically 
fit very well. In case of the classic SMA, the 
SK spectral data lack the predicted distortion. 
As a consequence, the region gets shifted to- 
wards smaller mixing, and even then the SK 
zenith angle spectrum considerably (about 3.5 
a) disfavors this fit. Due to the reduced mix- 
ing angle, the predicted SNO charged-current 
rate is more than 2a above the measurement. 
Consequently, small mixing is excluded at the 
3.5a level. 

Figure 4 shows a close up view of the LMA 
region. The combined fits described above 
treat the 8 B flux free parameter, which 
is only constrained by the solar neutrino ex- 
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Table 3 

Parameters for the best fit points. The probabilities given in the fourth row of numbers are based 
on the difference of % 2 (with respect to the minimum). The six rows below show the six independent 
parts of the fit: the Ax 2 from a fit to the shape of the SK zenith spectrum and the five interaction 
rates (deviation probabilities are given in units of Gaussian standard deviation a which contains only- 
experimental uncertainties). The last four rows show the values of the minimized fit parameters. The 
8 B and the hep fluxes are free, while the 8 B neutrino spectrum shift as well as the SK energy scale and 
resolution shifts are constrained within the systematic uncertainties. 



Solution 


Large Mixing 


Quasi- Vacuum 


Low Am 2 


Small Mixing 




Angle (LMA) 


(Quasi-VAC) 


(LOW) 


Angle (SMA) 


Am 2 


6.9xl(T 5 


6.68xl0~ 10 


7.2xl0~ 8 


6.6xl0~ 6 


tan 2 6» 


0.38 


1.5 


0.66 


0.0012 


X 2 (46 dof; p x2 [%]) 


43.5 (57.7) 


53.5 (20.9) 


52.5 (23.7) 


58.9 (9.6) 


A X 2 (2 dof;p Ax2 [%]) 


0.0 (100.0) 


10.0 (0.7) 


9.0 (1.1) 


15.4 (0.05) 


&XSK CPA x 2 H) 


3.0 ( 1.2a) 


1.0 ( 0.5a) 


5.6 ( 1.9a) 


15.1 ( 3.5a) 


Ga Rate [SNU] 


71.7 (-0.6<7) 


64.4 (-2.1a) 


64.9 (-2.0cr) 


78.3 (+0.7a) 


CI Rate [SNU] 


2.99 (+1.9cr) 


3.09 (+2.4(j) 


3.03 (+2.1(7) 


2.43 (-0.6(7) 


SK Rate [%SSM] 


46.1 (-0.47) 


45.0 (-1.2(j) 


45.5 (-0.8(7) 


45.7 (-0.6(7) 


SNO CC [%SSM] 


34.0 (-0.4cr) 


38.9 (+1.8(7) 


37.7 (+1.3a) 


40.0 (+2.3a) 


SNO NC [%SSM] 


106 (+0.4a) 


76 (-2.0a) 


87 (-l.la) 


89 (-l.Oa) 


<f>s B [10 6 /(cm 2 s)] 


5.33 (+0.3(7) 


3.83 (-1.5a) 


4.41 (-0.8a) 


4.48 (-0.7a) 


<f> hep [10 3 /(cm 2 s)] 


36 


9 


23 


12 


8 B Spectrum Shape 


-0.3(7 


-0.3(7 


+ 0.0(7 


+0.9(7 


SK E-scale/resol. 


-0.5ct/-0.1c7 


-0.5a/ 0.2a 


+0.0cr/-0.2cr 


+1.3a/-0.3(7 



periments. If the SSM prediction of this flux 
is used (dashed-dotted line), the quasi- VAC, 
LOW and small mixing solutions are further 
disfavored: the quasi- VAC solution is then ex- 
cluded at 99.4% C.L., LOW at 99.1% C.L. and 
small mixing at 99.96% C.L. If the SNO mea- 
surements are removed from the fit (dashed 
lines), LMA is still favored. Quasi-VAC solu- 
tions appear then at 95% C.L, LOW solutions 
at 97% C.L. and small mixing solutions at 
93% C.L. If the radio-chemical measurements 
(Gallium and Chlorine) are ignored (dotted 
lines, or see Figure 3 a), LOW solutions ap- 
pear at 92% C.L., quasi-VAC at 84% C.L. 
and small mixing solutions at 99.99% C.L. 



4 Conclusion 



The combined fit to the charged- current inter- 
action rates measured by several experiments 
result in many allowed regions of neutrino 
mixing and mass square difference: LMA, 
SMA, LOW, quasi-VAC and VAC solutions. 
The absence of spectral distortion and daily 
variations of the SK solar neutrino interac- 
tion rate strongly constrains these regions 
rejecting SMA, LOW and VAC solutions at 
95% C.L., while leaving only the higher mass 
square difference LMA and the quasi-VAC re- 
gion. When combined with the SK interaction 
rate and either the SSM prediction of the 8 B 
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flux or the SNO rates, large neutrino mixing 
is selected. When Gallium, Chlorine and SNO 
data are combined with SK data, only the 
higher mass square difference LMA solutions 
remain at 98.9% CL. Thus the combined re- 
sults of all solar neutrino experiments can be 
used to determine a unique region of oscil- 
lation parameters that explains the famous 
solar neutrino problem. 
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[%\ A detailed discussion of this uncertainty will 
be published elsewhere. 

|j| As we prepared this paper, SNO announced 
new results on its charged- and neutral- 
current rates, which assumes an undistorted 
8 B v e spectrum. In particular, the neutral- 
current rate can be strongly influenced by 
such spectral distortions. At present, the 
SNO measurement of this spectrum cannot 
be used to extract neutral-current and 
charged-current rates for distorted spectra, 
since the systematic uncertainties have not 
yet been given as a function of energy. As 
a consequence, we do not use the neutral- 
current rate in oscillation fits excluding SK 
data. In combined fits including the SK recoil 
electron spectrum, the distortions are small 
in the regions which are allowed by SK data. 
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